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THE BI-POISSON PROCESS: A QUADRATIC HARNESS 
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University of Cincinnati, Warsaw University of Technology 
and Warsaw University of Technology 

This paper is a continuation of our previous research on quadratic 
harnesses, that is, processes with linear regressions and quadratic 
conditional variances. Our main result is a construction of a Markov 
process from given orthogonal and martingale polynomials. The con- 
struction uses a two-parameter extension of the Al-Salam-Chihara 
polynomials and a relation between these polynomials for different 
values of parameters. 

1. Introduction. The family of stochastic processes with uncorrelated 
increments, linear regressions and quadratic conditional variances includes 
the Wiener, Poisson, Gamma and Pascal processes; for details and additional 
references, see [20], Theorem 1. It includes their free counterparts which are 
Markov processes whose bivariate distributions match the bivariate distri- 
butions of free Brownian motion and free Poisson process; for details and 
additional references, see [4], [10], Theorem 4.3, and [5], Proposition 3.4. It 
also includes the classical version of the q-Brownian motion, see [6], and the 
classical version of the q- Poisson process, see [1]. All of these examples share 
the property that the conditional variance with respect to the past sigma 
field is constant. 

In this paper we study the three-parameter family of stochastic processes 
with linear rather than constant conditional variances with respect to the 
past or the future sigma fields. The bi-Poisson process with parameters 
(r],8,q) is the square-integrable stochastic process (X t ) such that for all 
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< s < t < ii we have 

(1) nX t ) = 0, E(X t X s ) = mm{t, s}. 

(2) nX t \F a , u ] = —X s + —X u 

u — s u — s 

t 7- Iy , r i (« -t)(t-s) ( uX s -sX u 

Var I t f s , u = 1 + ?7 

u — qs V u — s 



(3) + 



■u — s 

(uX s -sX u )(X u -X s 



(1-9) 



(u — s) 2 
where 

.F s , u = a {X t : < t < s or f > u}. 

Thus in the terminology of [8] , the bi-Poisson process is a quadratic harness 
with parameters r = a = 0. 

We use the standard g-notation 

[n] q = l + q+--- + q n -\ 

[nW. = [l] q [2] q ---[n] q , 



[n-k\ q \[k] q V 

with the usual conventions [0] q = 0, [0] q \ = 1. 

In [8], Proposition 4.13, we give the following uniqueness result. 

Proposition 1.1. Fix — 1 <q< 1. If polynomials {p n (x;t)} given by 
xp n (x;t) =p n+ i(x;t) + (0 + tf7)[nLp n (a;;£) 

(4) 

+ 1 (1 + r}9 [n - 1] ,) [n] g p n _i (x ; t) , 

n > 0, with p_i = 0,po = 1; are orthogonal martingale polynomials for a 
Markov process (Xt), then (Xt) is a bi-Poisson process, that is, (2) and 
(3) hold, and \ + r]9> max{g, 0}. Moreover, such a Markov process (Xt) is 
determined uniquely. 

The goal of this paper is to construct the three-parameter family of bi- 
Poisson processes within the full range of the parameters. 

Theorem 1.2. For —l<q<l and 1 + r\9 > max{g, 0} there exists 
a Markov process (Xt) with Xq = and transition probabilities such that 
{p n (x;t)} are martingale and orthogonal polynomials for (Xt). 
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Combining Proposition 1.1 with Theorem 1.2 we get the following. 

Corollary 1.3. For every —1 < q < 1, 1 + r]6 > max{g,0}, there exists 
a unique bi-Poisson process (Xt) with parameters (rj,9,q). Moreover, (Xt) 
is Markov. 

Polynomials (4) are the so-called Al-Salam-Chihara polynomials and have 
been studied by many authors; see [13], page 377, [14], Section 3.8. A nat- 
ural plan of constructing the Markov process corresponding to a prescribed 
family of orthogonal polynomials relies on the Mehler-Poisson formula for 
the polynomials; such a plan was followed successfully in [7] where the poly- 
nomials were simple enough with rj = 9 = 0. However, the Mehler-Poisson 
formula for the more general polynomials (4) is not available. We therefore 
follow the approach from [10], where we define the Markov transition prob- 
abilities of (Xt) using auxiliary orthogonal polynomials {Qn}, which were 
guessed by a trial and error method. This approach was then tested on the 
bi-Poisson process with parameter q = in [12], where polynomials {p n } 
and also {Q n } satisfy a constant coefficients recurrence. The difficulty of 
the general bi-Poisson case lies in the fact that the corresponding orthogo- 
nal polynomials {Q n } satisfy a more complex three step recurrence which 
falls outside of the classes studied in the literature. The proof of the Markov 
property relies on nontrivial identities for connection coefficients between 
these orthogonal polynomials; the form of the identities deviates signifi- 
cantly from their predecessors [10], Lemma 3.1, [12], Proposition 2.2, and 
[9], Theorem 1; the latter identity expressed the Al-Salam-Chihara polyno- 
mials in terms of the continuous g-Hermite polynomials and provided the 
initial breakthrough that made our constructions feasible. 

Note added late. After the first version of this paper was written, we 
learned about the so-called "projection formulas" for the Askey-Wilson and 
related polynomials in papers by [17, 18, 19]. From the results in these papers 
one can deduce that there is a family of probability measures {P s j(x, dy)} 
such that the integral representation 



and for all x from the interval (15). One can further check that the Chapman- 
Kolmogorov equations hold for s < t from the interval (5), so measures 
{P St t(x, dy)} are the transition probabilities of the Markov process from 




holds for all t <u, all u from the (possibly empty) interval 
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Theorem 1.2. Thus one can deduce the existence of the process (Xt) for t in 

(5) . 

The proofs in [17, 18, 19] rely on infinite product identities for hyperge- 
ometric functions and do not cover the full range < t < oo. On the other 
hand, Theorem 1.2 does not generalize [17, 18, 19], as in our setting two of 
their parameters vanish. 

2. Orthogonal polynomials and the bi-Poisson process. We will deter- 
mine the transition probabilities of the general bi-Poisson process in implicit 
way as the orthogonality measures of the auxiliary six-parameter family of 
orthogonal polynomials {Q n (y;x,t,s) :n > —1} which satisfy the three step 
recurrence 

yQ n (y,x,t,s) 

(6) = Q n+1 (y;x,t,s) + A n {x,t, s)Q n (y; x,t, s) 

+ B n (x,t,s)Q n -x(y;x,t,s), 
with Q-i = 0, Qo = 1 and with the coefficients 

(7) A n (x, t, s) = q n x + [n] q (tr) + 9- [2} q q n ~ l srj) , 

(8) B n (x, t, s) = [n] q (t - sq™- 1 )^ + r^xq^ 1 + [n - l} qV (9 - s^- 1 )}, 

where n > 0. [For n = 0, the formulas should be interpreted as Ao(x, t, s) = x, 
Bo(x,t,s) = 0] 

In order for such a definition of transition probabilities to make sense, 
we will need to establish a number of properties of these polynomials. First, 
we need to know that the orthogonality measure exists for all the values of 
parameters we need. Second, we will need to show that the measures are 
consistent, that is, that the Chapman-Kolmogorov equations are satisfied. 
Finally, we need to verify that properties (1), (2) and (3) hold. 

For k = 0, 1, . . . , n, j = 0, 1, . . . , k let 

(9) In k J = sWq^-i)^, ~ k+J - 



[n-k] q \[n-k-l] q \[j] q \[k-j] q V 
with the conventions that: 

• 7n,fc,0 = 



, 0<k<n, 

• ln,n,j = for < j < n, while 7 nj n i0 = 1, 

• 7n,fcj = if j < 0, or k < 0, or j > k, or k > n. 

Define auxiliary polynomials 

k 

(10) bl n) (y;x,s) = ^2^ n ,k,jQk-j(y;x,0, s), k = 0,l,...,n. 

j=0 
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In particular, bn (y;x,s) = Q n (y;x,0,s), b^ \y;x,s) = 1 and b ™ (y;x,s) = 
if j < or j > n. 

We now state the main result of this section. 

Theorem 2.1. If polynomials {Q n } o,re defined by recurrence (6) with 
coefficients (7), (8) and {b k } are given by (10), then 

n 

(11) Q n (z;x,u,s) = Y^b^l k (y;x,s)Q k (z-,y,u,0) 

k=0 

holds for n = 1, 2, . . . . 

We remark that in fact a more general version of (11) holds. Let 
~ h .( t )-(-ny [n] q \[n-k + j-l] q \ ^ 

— k 

h { n) x ' *' s ) = Y InAjitfQk-j (y,x,t,s), k = 0, 1, . . . , n. 

3=0 

Then 

n 

Q n (z; x, u,s) = ^2 %l k (y; x, t, s)Q k (z; y, u, t). 

k=0 

The latter formula specializes to (11) when t = 0. For our purposes, (11) 
with y = would suffice. However, our argument relies on the three step 
recurrence with respect to variable y; this would become a mysterious trick 
if y were set to zero, so we will prove (11). 

Corollary 2.2. 

n 

(12) Q n (z]x,u,s) = ^b^} k (0;x,s)(p k (z;u) - p k (x]s)), n > 1. 

k=i 

Proof. Since Q n (z;0,u,0) =p n (z;u), clearly (11) with y = implies 

Q n (z;x,u, s) = Y,k=o^n-k(^ x ' s )Pk( z '^ u )- From (6) we see that Q n (x;x,s,s) = 

0, so applying the same formula again we get J2k=o ^n-k(^' x i s )Pk{ x ] s ) = 0. 
Subtracting these two sums and noting that po = 1 we get (12). □ 

The following proposition reduces Theorem 2.1 to the identities between 
{b^ k (y;x,s)}. 
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Proposition 2.3. Suppose a family of polynomials {Q n (y;x,t,s)} is 
given by recurrence (6) with arbitrary coefficients A n and B n such that 
Bo(x, t, s) = 0. If {b k (y; x, s) :n = 0, 1, . . . , k = 0, 1, . . . , n} satisfy equations 

bf +l \y;x,s) = b ( ^} 1 (y;x,s)(A n+1 - k (y,u,0) -A n (x,u,s)) 

+ b k %{y, x, s)B n+2 -k(y, u, o) 

(13) 

- b k-2 (y; x ' s ) B n{x, u, s) + (y; x, s), 

fc = 0,l,2,...,n + l, 

with the initial conditions b^ = 1 and b k = for k > 1 , £/ien (11) holds 
true. 

Proof. We proceed by induction on n. Obviously, (11) holds for n = 
and it is also easy to verify for n = 1. Suppose now that n is a positive 
integer such that (11) holds for all positive integers not greater than n. 

By (6) and inductive assumption 

Q n+ i(z;x,u,s) 

= zQ n (z; x, u, s) - A n (x, u, s)Q n (z; x, u, s) - B n (x, u, s)Q n -i(z; x, u, s) 

n 

= Z Y1 b n-k(v'i x > s )Qk(z; y, u, 0) 

k=0 

n 

- A n (x, u, s) b n-k(y'' x > s )Qk( z ; y, u, o) 

k=0 
n-1 

- B n (x,U,s) ^ b n-l-k(y^ X ^ S )Qk(z-,y,U,0)- 

k=0 

Again by (6), the right-hand side equals 

n 

XI h n-k^ x ' s ){Qk+i(z; y, u, 0) + A k (y, u, 0)Q k (z; y, u, 0) 

k=0 

+ B k (y, u, 0)Q k -i(z; y, u, 0)} 

n 

- A n (x, u, s) ^ b^}_ k {y- x, s)Q k (z; y, u, 0) 

fc=0 
n-1 

-B n (x,u,s) ^ b^} k (y;x, s)Q k (z-,y,u,0). 

k=0 
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Collecting the coefficients at Qj(z;y,u,0) we get 
Q n+1 (z;x,u,s) 

n+l 

= Qn-k+i(z;y,u,Q){bl n \y;x,s) 

k=0 

+ An-k+i (y, u , o) & i-i (y ' x > s ) 

+ Bn-k+2(y, u, 0)b^} 2 (y; x, s) 

- An(x,u, 8)6^(2/; x, s) 

- B n (x, u, s^Jp (y; x, s)}. 

Thus (11) follows from the assumed recurrences for &i n (recall the 
convention b^ 1 = when j < or j > n + 1). □ 

Proof of Theorem 2.1. By Proposition 2.3, in order to prove Theo- 
rem 2.1 it suffices to verify that functions {b^} as defined by formula (10) 
satisfy recurrence (13). This in turn consists of expanding both sides of (13) 
using (10), then applying (6) to Qj(y;x,0,s), and finally matching the co- 
efficients at Qj(y; x, 0, s). The last task is tedious so we postpone it to the 
Appendix. □ 

We now use Theorem 2.1 to prove our main result. The case q = is 
covered by our argument when fc is interpreted correctly as or 1; how- 
ever this case is much simpler and is available separately in [12]. The case 
q = —1 is elementary and is handled separately in Section 3.2. The case 
q = 1 requires additional modifications; since an elementary self-contained 
construction based on birth/death processes is available in [11], in Section 
3.1 we identify the transition probabilities and omit the rest of the proof. 

2.1. Proof of Theorem 1.2 for < \q\ < 1. If 07] = then (X t ) is a 
q-Brownian motion when both parameters are zero, a q-Poisson type process 
when rj = 0, 9 ^ 0, or a time inverse of a g-Poisson type process when 6 = 0, 
rj^O, see [10]. 

We will be therefore interested only in the case 0i] ^ 0. Passing to (—Xt) 
if necessary, without loss of generality we may assume T] > 0. Replacing our 
process by the process (\/v/\^\^t\e\/r])t>0i without loss of generality we may 
further assume that 6 = dbr/. 

The plan of the proof is similar to the proof for q = as given in [12]. 
We first verify that B n (x,s,t) > on the support of X s , so that we can de- 
fine the transition probabilities P s j(x,dy) as the orthogonality measures of 
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{Qn(y] x, t, s)}. We will then verify that these probability measures P s j(x, dy) 
satisfy the Chapman-Kolmogorov equation. Since p n (x;t) = Q n (x; 0, t, 0), 
the second property will imply that {p n (x;t)} are orthogonal martingale 
polynomials. 

Through the remainder of the proof we fix 

(14) \q\ < 1, 7] > 0, and such that 1 + r]0 > max{g,0}. 

2.1.1. Support of X t . Under assumption (14) the three step recurrence 
(4) has positive coefficient at p n -ii and the coefficients of the recurrence 
are bounded uniformly in n. Therefore, see [13], Theorems 2.5.4 and 2.5.5, 
polynomials {p n (x;t)} are orthogonal with respect to the unique compactly 
supported probability measure irt(dx). 

We determine 7r t by setting b = c = in [3], Theorem 2.5. To cast our 
recurrence (4) into the appropriate form, we need to re-parameterize our 
polynomials. With t > being fixed, let P n {x) = p n ((ax + P);t)/a n , where 

" = f^vV + l-9> P = T=f- Then (4) translates into 

xP n {x) = P„+i(x) + \(a + d)q n P n (x) adq n - l ){\ - ( ? ri )P ri _i(x), 

with 

i]0 1 9 + tr] 

ad = — , a + d — 



Tfl + l-q' y/i y/rjO + 1 - q' 

Solving this system of equations we get 

m&x{tr/, 0} 



d 



yfty/rje + l-q 
mm{trj, 0} 



y/ty/rjd + l-q' 

Thus {P n } are the Al-Salam-Chihara polynomials. 

Taking into account our linear change of variable ax + 0, from the general 
theory of Askey- Wilson polynomials [3], Theorem 2.5, see also [2], we see 
that the support of the absolutely continuous part of 7r t is the interval \J3 — 
a, (3 + a], that is, 

-0 + tJ]-2^i 1 + l-q + tr] + 2yV^ + 1 

(15) 



l-q 1-q 

for all t > 0. (The absolutely continuous part exists for all t > except when 
r]6 + 1 — q = 0, or when q = 1.) 

For < t < ri Q+ 1 _ q , the discrete part of irt is supported on the finite 
number of points 



(16) Xk = - ^Qq* + t > ' " -(tn + 0) 
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for k = 0, 1, . . . such that t(r)8 + 1 - q) < q 2k 6 2 . 

For t > 7?e ^j~ g , the discrete part of 7Tj is supported on the finite number 
of points 

(17) Xk = -J—( tm k + ^ + - (ty + 0) 



1 — q \ rjq K 

for A; = 0, 1, . . . such that tr] 2 q 2k > r]9 + 1 — q. 
These calculations lead to the following. 

Lemma 2.4. For fixed < t < u and B n given by (8) let 

(18) U t = f]lx:l[B n (x, U ,t)>0y 

n=l I k=l J 

Then irt(Ut) = 1. 

(Notice that set Ut does not depend on variable u.) 

Proof. We first verify that the absolutely continuous part of nt charges 
all of its mass to Ut- To see this, we show that it assigns no probability to 
any of the sets V n = {x : 1 + nq n ~ 1 x + [n — l] q r](9 — tnq n ~ l ) < 0}. Since B n is 
a linear function of x, to verify the latter, we only need to verify that none 
of the endpoints of the interval (15) is in V n for all n > 1, that is, 

1 + VQ + [n - l] q r](9 - tr/q n L )>0. 

We simplify this inequality to 

(rfi + l-q) + tn 2 q 2n - 2 > ±2 m n ~ 1 ViVr]0 + 1 - q, 

which holds true by the elementary inequality A 2 + B 2 > ±2AB. 

We now verify that the discrete part of 7r t does not charge the set XJ[ = 
U^Li V n . If t = then ir = 5 Q , so Tr (V n ) > implies 1 + r]6[n - l] q < 0, 
contradicting (14). Suppose t > and x^ is in the support of the discrete 
part and is given by (16). Then t(r]9 + 1 — q) < q 2k 9 2 and the condition 
£> n+ i(xfc,u,t) > is equivalent to 

(19) 1 + m n x k + [n] q r){9 - trjq n ) > 0. 

A calculation shows that this inequality is equivalent to 

1 - (i _ q + n e(\ - q n + k )) > 0. 



Notice that 1 — q + r)0(l — q n+k ) = (1 — q)(l + r]9[n + k] q ) > 0, as the parame- 
ters are chosen to ensure that the coefficient at the third term of recurrence 
(4) is nonnegative. So to show that (19) holds we only need to verify that 

q n ~ k tri 
1 - „ > 0. 



10 



W. BRYC, W. MATYSIAK AND J. WESOLOWSKI 



If q n k 9 < then both terms are positive and (19) holds. If q n k 9 > then 
t < q 2k 9 2 / (r}9 + I — q) implies 

1 - t^L > ! - ***-f£- = (1 - ,) 1 + + > 0, 
r]0 + 1 — q r)6 + 1 — q 

so (19) follows. 

Suppose now that Xk is in the support of the discrete part and is given 
by (17). Then the condition B n +i(xk,u,t) > is equivalent to 



1 — n n ~ k 

(20) =— ? {l-q + r,e-tr, 2 q n+k ) > 0. 

1 - q 

In particular, Bk+i(xk,u, t) = 0, so we only need to show that B n+ i(xk,u, t) > 

for < n < k. 

Assuming n < k we consider two cases. Suppose first that k — n > is 
even or q > 0. Then the first factor in (20) is negative and the second factor 
is (—tr] 2 q 2k + 1 — q + r]6) + rj 2 tq 2k (l — q n ~ k ) < as the sum of two negative 
terms. Suppose now that k — n > is odd and — 1 < q < 0. Then both factors 
in (20) are positive. □ 

2.1.2. Transition probabilities. This part of the proof follows closely [10] 
and [12]. For < s < t consider the six-parameter family of polynomials 
{Qn(y, x, t, s) :n > —1} in variable y which is defined by the three step re- 
currence (6) with the coefficients (7), (8). (In this notation, we suppress the 
dependence of Q n on the parameters rj, 9, q, which are fixed.) 

We define P s j(x,dy) as the (unique) probability measure which makes 
the polynomials {Q n {y\ x,t, s) :n £ N} orthogonal. This is possible whenever 

1 + r\9 > max{<7,0}, — 1 < q< 1, and x £ U s as defined by (18); the latter 
holds true for all x from the support of the probability measure it s (dy) = 
Po,s(0,dy), see Lemma 2.4. Since the coefficients of the three step recurrence 
(6) are bounded, it is well known that measures {P St t(x, dy)} have bounded 
support. 

The proof of Favard's theorem as given in [13], Section 2.5, shows that 
P Si t(x,dy) is a weak limit of the discrete measures supported on the zeros 
of polynomials {Q n (y;x,t,s)}; since these zeros depend continuously on x, 
functions x i— > P s j(x, A) are measurable. 

We now follow the argument from [10], Proposition 3.2, and verify that 
probability measures {P St t(x,dy)} are the transition probabilities of a Mar- 
kov process. 

We begin by noting that since / Q n (z; x, u, s)P S)U (x, dz) = for n > 1 and 
&o n ^(y;x,s) = 1, from (12) we recurrently get 



(21) J p n (z;u)P StU (x,dz) =p n (x), 



n > 1. 
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Proposition 2.5. IfO < s <t, l + ryfl > max{(7, 0} ; -1 < q < 1, then for 
a set of x of tt s -measure one, 

(22) P«,«(av)= / ?t,u(v, 

JUt 

Proof. We first verify the special case s = of (22), where we rename 
the variables for ease of use in the remaining proof: 

(23) n t (A)= f P S)t (z,A)n s (dz). 

J Us 

To see this, let v{A) = f v P Sj t(x,A)-K s (dx). Since P S) t(x, •) is a well-defined 
non- negative measure for all 1 £ to verify (23) we only need to show that 
the recurrence (4) defines polynomials {p n (x; t) : n > 0} which are orthogonal 
with respect to v. Since ir s (U s ) = 1, this follows from (12) in a manner similar 
to the general case below, so the argument is omitted. 

Let V Sjt = {x£ U s :P S)t (x,U t ) = 1}. From (23) we see that 7r s (V Stt ) = 1. 
Indeed, suppose ir s (V S! t) < 1 so that TV s (U s \V Sl t) > 0. Since P Sj t(x, Ut) < 1 on 
U s \V St t, we have 1 = ir t (U t ) = ir s (V S)t ) + J Us \ Vs t P St t(z,U t )7C s (dz) < TT s (V S)t ) + 
7T S (U S \ V Si t) = 1, a contradiction. 

Let 

u(A)= [ P tyU (y,A)P S)t (x,dy). 
JU t 

We will show that v{dz) = P S)U (x,dz) for all x G V Si t by verifying that the 
polynomials {Q n { z \ x, u , s)} are orthogonal with respect to v(dz). Since poly- 
nomials {Q n } satisfy the three step recurrence (6) with bounded coefficients, 
they determine a unique probability measure; since u(-) is a probability mea- 
sure, to verify that it coincides with P StU (x,dz) it suffices to show that for 
n > 1 measure v(dz) integrates Q n (z;x,u,s) to zero. 

We use the fact that f R Q n (z;y,u,t)Pt t u(y,dz) = for n > 1. From (12) 
and (21) we get 

Q n (z;x,u, s)u(dz) 
= / 5Z & i-fc(°; x ' s ) / (Pk(z;u) -Pk(x;s))P tu (y,dz)P st (x,dy) 

JUt JR 
n „ 

= J2 b n-k(°' x ^ s ) / (Pk(y,t) -p k {x;s))P S! t{x,dy) 
k=i JUt 

n . 

= J2 b n-k(°'> x ' s ) / (Pk(y,t) -Pk(x;s))P s j(x,dy) = 0, 
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as for x £ V s> t we have P S) t(a;,M \ Ut) = 0. □ 

The following lemma must be known but we could not find an appropriate 
reference. 

Lemma 2.6. If the family of probability measures {P Sj t(x, dy) :0 < s < t} 
satisfies (22) on the set of x's of Po, s (0, dx) -measure one, then there is 
a Ut-valued Markov process (Xt) with Xq = and transition probabilities 
P Sjt (x,dy). 

Conclusion of the proof of Theorem 1.2. By Proposition 2.5 and 
Lemma 2.6, there is a Markov process (Xt) such that Xq = with transition 
probabilities P St t(x,dy). 

Since p n (x;t) = Q n (x;0,t;0), the construction shows that {p n (x;t)} are 
orthogonal polynomials with respect to measure 7Vt(dx) = Po,t(0, dx). For a 
Markov process, the martingale polynomial property follows from (21). □ 

3. Some special cases. In this section we complete the proof of Theorem 
1.2 by analyzing two previously omitted extreme values of q. 

3.1. Bi-Poisson process with q = l. In this section we read out the tran- 
sition probabilities determined by (6) in the case q = l; then the remaining 
parameters satisfy nO > 0. This case is not covered by the construction in 
Section 2.1 as the coefficients of recurrence (4) are unbounded. After we 
identify the transition probabilities from (6), we will omit the construction, 
as an even more elementary construction based on birth/death processes is 
presented in [11]. 

Proposition 3.1. Suppose n,9 > and q = l. If (Y t ) is a Markov pro- 
cess with transition probabilities such that polynomials defined by the three 
step recurrence (6) with coefficients (7), (8) are orthogonal, then 



Y t 



(6- V t)Zt-t, 0<t<9/ri, 

Z e/ --, t = 6/rj, 

(r)t-6)Z t -~, t>e/ V , 
V 



where (Zt) is Markov with the transition probabilities determined as follows. 
Z t -Z s \Z s =NB^ + Z s ,^-^j, O<s<t<0/ V , 

Z e/v \Z s = Gam(^- + Z s ,e-7 1 s^, 0<s<6/ri, 
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Table 1 

Distributions of the bi-Poisson process with q = 1 



Name Parameters Distribution E(e uZ ) Notation 

Poisson A>0 e~ x \ k /k\, k = 0, 1, 2, . . . exp(A(e" - 1)) Poiss(A) 

Gamma r > 0, a > f(z) = ^n^ e~ x/a (l-au)~ r Gam(r,cr) 

V j_ — p) , k — 



Negative r > 0, 0<p<l T^p r (l - p) k , k = 0, 1, . . . a-a-PWV NB(r,p) 
Binomial 

Binomial n > 0, 0<p<l (£)p fe (l -p)"~ fe , /c = 0, 1, . . . , n. (l-p + pe 11 )' 1 Bin(n,p) 



Z^/^Poissf-^), i>0/r?, 



\Tft 

Z t \Z s = Bm^Z s ,^^j, 9/ V <s<t. 

Here we use the fairly standard notation for the distributions summarized 
in Table 1. 

The following result can be traced back to [15]; we state it in the para- 
metrization used in [10]. 

Lemma 3.2. Fix parameters t > 0, 9, f G M. Suppose polynomials {p n (y) '■ 
n > 0} are given by the three step recurrence 

(24) ypn{y) = p n +i(y) + 0np n (y) + (t + ?(n- l))np n _i(y) 

/or n = 1,2, . . . uiii/i p_i = and po = 1. 

Let Y be a random variable such that {p n (Y)} are orthogonal. 

(i) // f = 0, 9 + 0, i/ien 

y = #z - t/0, 

where Z is Poisson with parameter X = t/9 2 . 

(ii) If9 2 = 4f>0 then 

Y = sgn(9)Z - 2t/9, 

where Z is Gamma with shape parameter r = t/f and scale parameter a = 
|0|/2. 

(iii) // fl 2 > 4f > then 



'6 2 - Art 

Y = sgn(0W# 2 - AfZ - sgn(#) — v — 1, 

2f 

where Z is negative binomial with parameters r = t/f, p — ~ 4 — 
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(iv) If f < and n = — t/f is a natural number then 



Y = \/9 2 - 4fZ + 7^r(V# 2 - 4f 
2r 



where Z is binomial Bin(n, p) with p = i(l 



Note. If (24) holds with i = then Y = 0. 

Proof. The type of the law of Y can be read out from [16], Table 1 
and formula (8.3). Perhaps the simplest way to identify the parameters of 
the basic law of Z as well as the linear transformation Y = a(Z — E(Z)), is 
to compare the first four moments. We have 

E(F)=0, E(Y 2 ) = t, K(Y 3 ) = 9t, E(Y 4 ) = 3P + (9 2 + 2f)t. 

Solving the system of equations a J K(Z — KZ) J = E(Y J ), j = 2,3,4, we get 
the results as stated. 

Except for (iv), the type of the distribution can also be read out from [10], 
Theorem 4.2. However caution is needed as some of the explicit functions 
are typeset incorrectly in the published version (the formulas are corrected 
in the arXiv preprint). □ 

Proof of Proposition 3.1. Fixed < s < t. Recurrence (6) for poly- 
nomials Q n {z + x; x, t, s) is equivalent to (24) with the following identification 
of parameters: 

= trj + 6 - 2srj, f = r,(6 - srj)(t - s), t = (1 + rjx)(t - s). 

A calculation shows that 6 2 - 4f = (t-q - 9) 2 . Thus if Y s > — l/rj then t > 
and one of the cases listed in Lemma 3.2 must hold. 

We now use (6) with x = s = to verify that for t > we have Y t > — l/rj. 
For t^6/r] by Lemma 3.2 we have 

min{^, rjt} 



Y t = \6-Tft\Z t 



6r] 



where Z t ~ NB(^, m ^ { g^ } ). For t = 6/rj we have Y e/r) = Z - |, where Z ~ 
Gamma(^,^). Therefore, Y s > —1/r/ and the inequality is strict except when 
s > 9/rj. 

We now use Lemma 3.2 to identify the conditional distribution of Yt given 
Y 8 . 

(i) Suppose < s < t < 9/rj. Then 6 > 9 - trj > and f > 0. 
Y t - Y S \Y S 2(9- t V )Z - (t - s)\±^, 

9 — ST] 
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where Z\Y S is negative binomial with parameters r = Y g^(^ and p = gZ^l - 
Writing Y s = (6 — r]s)Z s — s/6 we get r = ^ + Z s and 

Y t - Y S \Y S = {6- t V )Z - (i - s)/9 - (t - s)Z 8 r } . 
Letting Z = Z t — Z s we get 

Y t \Y s = {6-tr,)Z t -t/6. 



(ii) Suppose < s < t = 9/r}. Then 6 2 - 4f = and f > 0. Thus Y t 

^Z 

Y s + l/7 



Y S \Y S = Zt — Y s — 1/t], so Yt\Y s = Z% — l/i], where Zt\Y s is Gamma with 



Xz+lhL a nd a = 9- srj. Writing Y s = {9- rjs)Z s - s/6 we get r = 4 + Z s . 



(iii) Suppose s = 6/r]< t. Then f = 0, 6 = trj - 6 > 0. Thus Y t - Y S \Y S = 
(t V - 6)Z t - (1/r? + Y s ), that is, Y t \Y s = (tr, - 6)Z t - 1/ V , where Z t \Y s is 
Poiss(A) withA = ^ = ^. 

(iv) Suppose 6/rj < s <t. Then f < 0. Therefore, there must be an integer 
n = Z s such that t = — fn, that is, Y s = (srj — 6)Z S — \/r\. Then Yt — Y S \Y S = 
(trj — 6)Z t — (Y s + 1/r/) = (trj — 6)Z t — {srj — 6)Z S , where Z t \Z s is binomial 
with p = j//E§ an d n = Z s . Thus 

Y t \Z s = (t v -6)Z t -l/ V . 

(v) Suppose < s < 6/rj < t. Then 6 = (t - s)rj + 9 - srj > and f > 0. 
Therefore Yt — Y S \Y S = (trj — 6)Z — 1+vYs ; where Z\Y S is negative binomial 
withr = 4±^andp=^y.Thus' ? 

Y t \Y s = (trj-6)Z-l/n. 
Writing Y s = (6 — rjs)Z s — s/9, Z = Z t , we get r = -K + Z s and 



Y t \Y s ^(tr 1 -6)Z t -l/ V . 



This shows that (Z t ) is a Markov process with the transition probabilities 
as stated. □ 



3.2. Bi-poisson process with q = — 1. In this section we read out the tran- 
sition probabilities determined by (6) in the case q = — 1; then the remaining 
parameters satisfy 1 + r\6 > 0. This case is not covered by the construction 
in Section 2.1, as the reasoning based on [3] fails. 

We first observe that the polynomials p n {x\ t) = on the support of irt for 
n > 1, t > 0. Indeed, since [2] 9 = 0, from (4) we see that p n (x;t) has ^(a^) 
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as a factor for n > 2. Thus the univariate distributions of (Xt) are supported 
at two points a±(t) which are the roots of p2(x;t). The roots are 



o±(t) 



tr, + 9 ± J At + (tr? + Of 



2 

From the requirement that pi(x;t) = x is orthogonal to P2 (i.e., that the 
mean of Xt is zero), we compute the corresponding probabilities 

p ± (t) :=Pr(X t = a±(t)) = ±T 



2 ^ 2^1J+W+W' 

A calculation confirms that K(X^) = a+p+ + a?Lp- = t. 

Similarly, see (6), the transition probabilities that make polynomials {Q n } 
orthogonal are supported on the roots of 

Q2(y,x,t,s) = y 2 + y({srj - x)(l + q) - tr] - 9) + s - t + qx 2 + x6 - qsxr/. 

A calculation shows that the supports of the one-dimensional distributions 
are consistent: 

Q 2 (a±(t);a±(s),t,s) = 0. 

We now determine the transition probabilities 

p±±{s,t) = Pr(X t = a±(t)\X s = o±(s)) 

from the requirement that ¥,(Xt\X s ) = X s . A calculation gives the transition 
matrix 



P,,t 



p—(s,t) p-+(s,t) 
p+-{s,t) p++{s,t) 



with 



p++(s,t) 
p-+(s,t) 
p+-(s,t) 
p—{s,t) 



sri-t7] + J4s + (srj + 9) 2 + Ju + {tr, + 0) 2 



2JAt + {tr, + ef 



sri-tri-J4s+ {sr, + 9) 2 + J4t+ {tr, + 9) 2 



2JAt+{tr, + 9f 



-sr] + tr,-J4s + (sr, + 9) 2 + <J4t+ (tr) + 9) 2 



2JAt + {tr, + 9) 2 



-sr, + tr, + J As + (sr, + 9) 2 + J4t + (tr, + 9) 2 



2J4t + (tr, + 9Y 
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Using these formulas, a calculation verifies the Chapman-Kolmogorov equa- 
tion 

Ps,tPt,u = Ps,u- 

Thus the two-valued Markov process X t G {a + {t), «_(£)} is well defined, and 
has mean zero and variance t. 

To verify (2) and (3), we compute the conditional probabilities: for a, f3, 

7 = ±, 

Pr(Aj = ap{t)\X s = a a (s),X u = a 7 (u)) = — . 

A long but straightforward calculation now verifies (2) and (3). 

Prom E[X t |X s ] =X S , we deduce E(X S X S ) = s for s < t so (1) holds. 

APPENDIX: COMPUTATIONS NEEDED FOR THE PROOF 
OF THEOREM 2.1 

Our aim is to show (13) for defined in (10) and (9). 

We first notice 6q = r y n ,o,oQo(y] x, 0, s) = 7 n ,o,o = 1, so (13) holds for k = 

0. 

It remains therefore to verify (13) for k > 1. We note that (13) takes a 
slightly different form for the boundary values k = 1, n, n + 1 however due 
to our convention these cases are covered by (13). 

Since the right-hand side of (13) is an affine function of u we split it into 
two equalities, which will be proved separately: the first for the coefficient 
of u 

[n + 2 - k] g {l + m n+l ~ k y + [n + 1 - k] qV e}b^} 2 (y; x, s) 

(25) =r ] [k-l] qq n+1 - k bi l } 1 (y;x,s) 

+ [n] q {i + m n ~ x x + [n - i] gV (0 - m n - l s)}btz 2 l) (y; x, a), 

and the second for the intercept 

b^ +1 \y;x,s) 

= {q n+i-k y+[n+1 _ k]qd 

(26) - q n x - [n] q (e - 77(1 + q)q n - 1 s)}bfl 1 {y- x, s) 

+ [n] q q n - x s{l + m n ~ l x + [n- l] q rj(6 - w™" 1 *)}^ fa x > «) 
+ b^\y;x,s). 

First step. We start with proving (25). We expand 6^!_\(y; x, s), b^}_ 2 {y, x, s) 
and b^_2\y;x,s) in (25) using (10) and then after applying the three terms 



18 W. BRYC, W. MATYSIAK AND J. WESOLOWSKI 

recurrence formula (6) we arrive at the following equality which is equivalent 
to (25): 

fc-l 

(27) c n ,k,iQj(y^As) = o, 

j=0 

where 

C n ,k,j = -r}[k - l] g q n+1 ~ k "fn,k-i,k-i-j 

+ 7][n + 2 - k] q q n+1 ~ k {'y n ^2,k-i-j + Aj(x, 0, s)j njk -2,k-2-j 

+ B j+1 (x, 0, s)7 ni fc_ 2 ,fc-3-j} 

+ [n + 2 - fc],{l + [n + 1 - A;] g r/6 , }7 nifc _ 2 ,fc-2-i 

- [n] q {l + T/q n ~^x + [n- l] q r)(0 - m™' 1 s)}ln~i,k-2,k~2-j- 

We will show that C n k j = 0. To this end we note that due to (7) and (8) 
the expression for C n j~j can be written as 

C n ,k,j = r]0[n + 2 - k] q a + q n+1 - k+ ^x[n + 2 - k] q b + q n+1 - k V c + d, 

where 

a = [j]qq n+1 ~ k ln,k-2,k-2-j ~ [j]q[j + l] q q j T}S^ k -2,k-Z-j 

[ n \q[ n ~ 1]<7 

+ [n + 1 — K „7 n k -2,k-2-j — 7 — r ri-7n-l,fc-2,fc-2-ij 

' [ra + 2 — fcj 

b = 7„ ;fe _ 2 ,fe-2-j - [j + l]< ? g J 'f/S7n,fc-2,fc-3-i 

Wo k-2-i 

-Q 7n-l,fc-2,fc-2-j) 



[n + 2 - jfe] 



c = — [fc — l] 7 nifc _i jfc _i_j + [n] g [n - l] 9 g fc 2 ??s7n-i,fc-2,fc-2 



+ [n + 2 - fc]g{7n,fc-2,fc-l-j - (1 + q)q 3 1 VSln,k-2,k-2 



-J 



+ bljj + l] g 2j (^) 2 7n,fc-2,fc-3-j} 



d = -[n + 2 - k] q \j + l],9 n+1 k+j rjsj n ^2,k-3-j + [n + 2- k] q <y njk - 2 ,k-2- 



3 



]g7n-l,fc-2,fc-2-j- 



We will show that a = b = c = d = 0. 
a = 0: From (9) we get 



[fc-2-i] 



(28) [j + l] g g : '??S7n,A : -2,fc-3-i = \ : TT-ln,k-2,k-2-j 

[n-l-j\ q 
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and 

(0Q , r i [n + 2-k] q [n + l-k] q 

[ZU) [n 9 7n-l ,k-2 ,k-2-j = F '-. Tj 7n ,k-2 ,k-2-j- 

[n-l-j\ q 

Consequently 

a = [ " + 2 - fe] f- fc - 2 - fc - 2 - J {[;W +1 - fc ([« - 1 - i] 9 - [k - 2 - j] ff ) 

l n 1 3\q 

-[n + l-k] q ([n-l] q -[n-l-j] q )}. 

Using twice the g-subtraction formula 

(30) [m] q -[l] q = q l [m-l] q , l<m, 

to [n — 1 — j] g — [fc — 2 — and to [n — l] q — [n — 1 — j] q we get a = 0. 
b = 0: To calculate b we first use (28) and (29), getting 



Now we use (30) to [n — 1 — j] g — [fc — 2 — j] q which gives b = 0. 

c = 0: For c we need two more formulas: equation (28) with j — 1 gives 

/oi\ [ n -j\q\j]q j-1 

(31) 7n,k-2,k-l-j = T, : TrQ ? ?«7n,fc-2,fc-2-j, 

and from (9) we get 

/oo\ [n + 2-fc] 9 [ra + l-fc] g fe _ 2 

(32) 7 ni fc_i ! fc_i_j = j- : — q rjsj n> k-2,k-2-j- 

Inserting (28), (31) and (32) into the definition of c we get 
_ ysq n - k+1 [n + 2 - fc]q7 n , fc -2,fc-2-j 



n 



1-jUk-l-j) 



q 



x{-[k- l] q [n + l-k] q [n- l-figq*- 1 -* 

+ [n ~ j]q[j]q[n ~ 1 ~j]q 

~ [j]q[ n ~ 1 - j]q[k ~ 1 - j]q ~ [j]q[n ~ 1 - j]q[k - 1 - j] q q 

+ \j]q[k-2-j] q [k-l-j] q q 

+ [n- l] q [n + l-k] q [k- l-j] q q k - l - j }. 

Using (30) for [n — j] q — [k — 1 — j] q we combine the second and the third 
expression in the brackets above, and applying (30) to [n— l — j] q — [k — 2 — j] q 
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we combine the fourth and the fifth expression in the brackets. Thus the 
expression in the bracket takes the form 

q*- 1 '^- [k - l] q [n + 1 - k] q [n - 1 - j] q + [j] q [n - 1 - j] q [n+ 1 - k] q 

- [j] q [k - 1 - j] q [n + 1 - k] q + [n-l] q [n + l- k] q [k - 1 - j] q }. 

Now we use again (30) to [k — l] q — [j] q to combine the first two expressions, 
and to [n — l] q — [j] q for the last two expression to conclude that c = 0. 
d = 0: First apply (28) and (29) to get 

d = [U + 2 = k ^:^ {[n - 1 - j] q - [n + 1 - k] q - q^~ k [k - 2 - j] q} . 
[n-l-j\ q 

Then by (30) applied to [n — 1 — j] q — [n + 1 — k] q we see that d = 0. 

Second step. Next, we prove (26). We expand b^ +1 \y;x,s), &i_ 1 (y;a;,s), 

6^ 2 1 ^(y;x,s) and b^\y]x,s) in (26) using (10) and then after applying the 
three terms recurrence formula (6) we arrive at the following equality which 
is equivalent to (26): 

k-l 

(33) D n,kjQj (y; x ' °> s ) = °> 

where 

Ai,fc,j = 7n+l,fe, fc-j — Q n+l k ln,k-l,k-j ~ Q n+l k Aj(x, 0, s)^j n ^-l t k-l-j 

- q n+l - k B j+1 (x, 0, s)-y n ,k-i,k-2-j 

- {[n + 1 - k] q 9 - q n x - [n} q {6 - (1 + g)g n ~ V)bnJ^i,fc-i-i 

- [n]^™" 1 ^ + J)^" 1 ! + [n - l]^ - V)}7n-i,fc-2,fc-2-i 

We will show that -D n ,fcj = 0. To this end we note that due to (7) and (8) 
the expression for D n ,k,j can be written as 

D n ,k,j = q n+1 ~ k+j xA + q n+1 - k 6B + q n+1 ~ k+ i sC + D, 

where 

A = -7n,fc-l,fc-l-j + [? + l]g<rV7n,fc-l,fc-2-j 

+ g fe_1_i 7n,ifc-l,fc-l-j - Ng9"" 3+fc " j ??S7 n _i ife _2,fc-2-i, 
B = -[j]qln,k-l,k-l-j + [j + l]g[j]gg J 7?S7 nifc _ ljjfc _ 2 -j 

+ [A; - l] g 7 n)fc _i i jfc_i_j - [n] g [n - l] <? 9 fc_2 r?s7 n _i ifc _2,fc-2-j, 
c = [j + l] ? 7„,fc-i,fc-2-j - [n] (? g fc " 2 "- J 7 n _i ifc _2,fc-2-i, 
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D = 7 n+ i, fcifc _ i - q n+1 ~ k Jn,k-l,k-j + [j'Wl + q)q n ~ k+j r]S-fn,k-l,k-l-j 

~ [j + l] q \j] q q^ 1 - k+2i (vs) 2 1n,k-l,k-2r- j 

- [n] q (l + q)q n ~ 1 r}S'y n ^-i,k-l-j 

+ [n] q [n - l]qq 2n ~ 2 (r]s) 2 "/ n - ltk -2,k-2-j ~ 7n,k,k-j- 

We will show that A = B = C = D = 0. 

A = 0: Equation (31) with k replaced by k + 1 gives 

[k — I — j] 

(34) [j + l^TjSJn fc _i ,k-2-j = 7 ; ^Tn.fc-l.fc-i-j- 

Moreover by (9) we have 

[k — 1 — j] 

(35) [n] q q k ~ 2 r]sj n - 1<k -2,k-2-j = 7 ; ^Tn.fc-i.fc-l-r 

Consequently, by (34) and (35), we get 

A = { - [n - 1 - j] g + [k-l- j] q 

+ q k ~ l -l [n-1- j] q - q^ [k - 1 - j],}. 

Applying now (30) to the first two terms in the brackets we get 

-([n-1 - j] q — [k — 1— j] q ) = -q k - l ~i[n - k] q . 

Combining this with the third term we get by (30) again 

g*-W( [7l - 1 - j] q - [n - k] q ) = q n ~^[k - 1 - j] q . 

Thus together with the fourth term we arrive at A = 0. 
B = 0: By (34) and (35) we get 

B = %r = T3^ {-ttU" - 1 - j\ + b1# - 1 - j], 

+ [fc - l] q [n - 1 - j], — [n — l] q [k — l- j] q }. 

Using (30) we have [n — 1 — j] q — [k — 1 — j] q = g fc_1_:, [n — k] q , which is 
applied to the first two terms in curly braces above and also to the third 
term. Thus the expression in curly braces takes the form 

-\j] q [n - k] q q k ~ x ~i + [k — l] q [n - k]^' 1 ^ 

+ [ k - l] q [k - 1 - j] q -[ n -l] q [k-l- j] q . 

Now combining the first two terms by (30) we get [n — k] q [k — 1 — i] g g fc_1 
and combining the last two terms, again by (30) we obtain — [k — 1 — j] q [n — 
k] q q k ~ l , and thus B = 0. 
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C = 0: By (9) it follows that 

[n] q q k ~ 2 -f n ~l,k-2,k~2-j = [j + l\ q q^ n ,k-l,k-2-j 

which implies immediately C = 0. 
D = 0: By (9), see (35), 

foa\ [ n + l ]q[ n -j]q k-i 

[ fc ~~ Jlq 

From (28) and (29) with k replaced by k + 1, we get, respectively 

[n — i]q[j}q j_i 

(37) 7„ ifc _i ifc _ i = ^ _ ^ — q J ??s7 nifc -i,fc-i-j, 

(38) W ^ = [n + 1 ~ k] f~ k]q q k - 1 r,s ln , k _ 1>k _ 1 _ j . 

[ fe ~ Jlq 

Plugging (34)-(38) to the definition of D we can express each term through 
-fn jk -i,k-i-j and thus 

-[n-j]q[j]q[n-l-j]qq n - k+ ^ 

+ \j]q[k-j]q[n-l-3\q{l + q)q n - k+] 

-\j] q [k-j] q [k-i-jU n+l - k+i 

-[n} q [n-l-3] q [k-j] q {l + q)q n - 1 
+ [n-l] q [k-j] q [k-l-j] q q 2n - k 
-[n + l-k) q [n- k] q [n - 1 - j}^' 1 }. 

Now we combine the fourth and the sixth term in the brackets using (30) 
three times: 

[k ~ J]qq n+1 ~ k+j {[n - l] q [k - 1 -JU"- 1 -* - [k - 1 - j] q [j] q } 

= [k-j] g q n+1 - k+j {m q + q j i n - 1 - jUk - i-JU"- 1 -* 

-{[n-l-j] q -q k - l -i[n-k] q )\j] q ) 
= [k- j] q {[n -l-j] q [ k -l- j] q q 2n +i- k - [n - 1 - j] q [j] q q n+1 ' k+ ^ 

+ \jy n ([n-k} q + q n - k [k-l-j] q )} 
= [k- j] q [n - 1 - j) q ([k - 1 - j] q q 2n - k+3 + [j] q (q n ~ q n+1 ~ k+J )). 
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After inserting this in the last expression for D we note that the factor 
[n — 1 — j)q appears in every summand in the brackets. Canceling it we get 

D = Wtf- 1 *- 1 -* {[n + lUn - j] q q^ -[n- j] q \j] q q^ 
l K ~ Jig 

+ mk-3] q q n - k+j + [k-j] q [k - l-j] q q 2n - k+i 
+ [k-j] q [jU n 

-[n] q [k-j] q q n - l -[n\ q [k-j] q q n 

- [n + l-k^ln-k]^- 1 }. 
Adding the second and the third term in the brackets above we get by (30) 

-[j] q q n - k+j ([n- j] q - [k-j] q ) = -[j] q [n - k] q q n 
and adding to this sum the eighth term we obtain 

-[n- k\ q q k -\[n + l — k] q + [j] q q n+1 ~ k ) = -[n- k] q [n + 1 - k + j] q q k ~ 1 . 
Now we combine this sum with the first term: 

q^Hn + l] q [n - j] q - [n - k] q [n + 1 - k + j] q } 

= q k - l {([n + l — k+j] q + q n+l - k+ i[k - j] q )[n - j] q 

-(["- j\ q ~ Q n ~ k [k ~ j] q )[n + l + k- j] q } 

= [k — j] q ([n - j] q q n+ i + [n + 1 - k + j] q q n ~ 1 ). 

Thus the above expression is the sum of the terms: first, second, third and 
eighth. Adding the remaining terms and taking common [k — j] q we get 

D = r}sj n ^-i,k-\-j q 17 " 1 ^- j] q q j+1 + [n + 1 - k + j] q 

+ [k-l- j] q q n+1 - k+j + [j] q q - [n] q - [n] q q}. 

Thus (30) implies that the sum of the first, fourth and sixth term is zero as 
well as the sum of the second third and fifth terms. Consequently D = 0. 
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